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Longitudinal spin diffusion of two pseudo-spin domains is studied in a trapped sample 

above quantum degeneracy, and the effect of coherence in the domain wall on the dynamics of the 
system is investigated. Coherence in a domain wall leads to transverse-spin-mediated longitudinal 
spin diffusion that is slower than classical predictions, as well as altering the domains’ oscillation 
frequency. The system also shows an instability in the longitudinal spin dynamics as longitudinal 
and transverse spin components couple, and a conversion of longitudinal spin to transverse spin is 
observed, resulting in an increase in the total amount of coherence in the system. 
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The development of ultracold trapped atomic systems 
of fermions and bosons has opened up new paths of study 
for transport phenomena. In particular, theoretical and 
experimental studies of spin transport have been of recent 
interest in a number of systems including strongly inter¬ 
acting degenerate Fermi gases near unitary [iMsI, dilute 
noncondensed bosonic gases close to quantum degeneracy 
Hi, and ultracold gases in optical lattices 0 . 

Quantum degeneracy effects can dramatically mod¬ 
ify spin transport in spin-polarized systems, for exam¬ 
ple leading to phenomena such as the Leggett-Rice effect 
[8[. Though spin is a fundamentally quantum property, 
classical Boltzmann theory was expected to describe spin 
transport in dilute nondegenerate gases accurately, since 
quantum degeneracy effects are typically insignificant. 
Surprisingly in 1982 Bashkin 0 and Lhuillier and Laloe 
0 independently predicted that transport in quantum 
gases can lead to macroscopic collective behavior even 
in nondegenerate systems with spin-independent interac¬ 
tions. The driving process for this collective behavior is 
the identical spin rotation effect (ISRE), in which sym- 
metrization of the two-body wave function in binary col¬ 
lisions between identical particles leads to coherent ex¬ 
change interactions when the thermal de Broglie wave¬ 
length is larger than the two-body interaction length 
[l0|. Examples of such collective behavior include spin 
waves 0, spin self-rephasing HH , and non-conservation 
of transverse spin [l2|. 


Diffusion is also strongly modified by quantum effects, 
leading to nonlinear and anisotropic spin diffusion 13, 1^ 
and universal spin dynamics [l|. Near unitarity, exper¬ 
iments on two spin domains in a strongly interacting 
Eermi gas showed reversal of spin currents 0 , result¬ 
ing in quantum-limited spin diffusion. Also in unitary 
Eermi gases, further quantum modifications to diffusive 
timescales were observed Hi- In 0, patterned spin 
textures were used to study diffusion in various geome¬ 
tries, showing diffusive behavior in one dimension (ID), 
but large departures from classical transport in 2D. 


In the nondegenerate case of a ID two-domain spin 
structure studied here, quantum modifications to dif¬ 
fusion appear as a decrease in the oscillation rate of 
spin domains in a harmonic potential and an increase 
in the longitudinal spin diffusion time, through coupling 
to the transverse spin. This effect shows sensitivity to 
the degree of coherence in the domain wall between spin 
domains, which highlights the quantum mechanical na¬ 
ture of the dynamics more clearly. In this Letter we 
demonstrate the alteration of longitudinal spin diffusion 
timescales depending on the initial amount of the coher¬ 
ence in the domain wall and also report observation of a 
longitudinal spin instability, which leads to an increase 
in the total amount of coherence in the system. 

The mechanism for coherent atom-atom scattering is 
the ISRE, where exchange scattering between indistin¬ 
guishable particles leads to precession of each atom’s spin 
about their combined spin. The theoretical description 
of spin dynamics in nondegenerate quantum sy^ems was 
developed using semiclassical kinetic theory [l3|, llbMlQj . 
We adopt the notation of 0 and describe the time evo¬ 
lution of the spin density distribution, a{p^z,t), with a 
ID quantum Boltzmann equation that includes the ISRE 
via a spin-torque term: 

dta{p,z,t) + doa{p,z,t) - x a{p,z,t) = (1) 

where and Q = (C/diffM|| + gS)/h. 

The mean-field interaction strength for two scattering 
particles with 5-wave scattering length a and mass m is 
g = ^ , and ffy is the longitudinal unit vector on the 

Bloch sphere. The experimentally observable quantity is 
the spin S{z,t) = f dpa{p, z,t)l27rh. S contains the lon¬ 
gitudinal spin S\\ and the transverse spin S_\_, with mag¬ 
nitude S± (i.e. the spin coherence). lAxt and I/diff are the 
trapping potential and differential potential experienced 
between the two states (due to mean-field and Zeeman 
shifts). Lastly, the right-hand side of Eq. [T] represents 
damping due to elastic collisions 0 . 
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We create a two-domain structure of “up-down” 
pseduo-spin states in a nondegenerate gas of atoms 
trapped in an external harmonic potential and observe 
the relaxation to equilibrium. The trap is axisymmet- 
ric, with trapping frequencies of 27r x (6.7,260,260) Hz. 
Due to rapid averaging from the high radial frequency, 
dynamics can be treated effectively as one-dimensional. 
The pseudo-spin doublet consists of two magnetically 
trapped hyperfine ground states of ^^Rb (|1) = |1,—1) 
and 12) = |2,1) ), c oupled via a two-photon microwave 
transition. (See (2^, [2l| for details.) 

We prepare spin domains with optical patterning and 
microwave pulse techniques (see Fig.[T]) [20|,l2l|. Magnet¬ 
ically trapped ^^Rb is evaporatively cooled to just above 
quantum degeneracy (temperature T ^ 650 nK and peak 
density n ^ 2.6 x 10^^ cm“^). Initially in state |1), a 
microwave 7r-pulse transfers atoms to |2), except where 
illuminated with an off-resonant, partially masked laser 
beam that creates an optical step potential on top of 
the trapping potential. The masked laser ac Stark shifts 
atomic energy levels out of resonance by ^ 50f^R. This 
procedure creates domains of spin up and down sepa¬ 
rated by a helical domain wall in which the spin vector 
is coherently twisted, while remaining fully polarized. 



FIG. 1. (a) Schematic of the domain preparation, showing 
the patterning laser and accompanying spin evolution in the 
Bloch sphere representation (longitudinal axis and trans¬ 
verse equatorial plane). Applying a 7r-pulse to atoms initially 
in |1) while illuminating them with the masked laser trans¬ 
fers atoms in the unilluminated half (left) of the distribution 
to 12) while leaving atoms in the illuminated half unchanged 
(right), (b) Composite image of atoms in |1) (purple) and |2) 
(green) following initialization. 

The longitudinal spin projection is measured by di¬ 
rectly measuring the population in states |1) and |2) 
using absorption imaging. Each image is divided into 
equally spaced axial bins, and the local population dif- 


ference is given by S'||(2:) = with = 

'^Ni{z), {i = 1,2) summed over all axial bins. Ramsey- 

Z 

type experiments measure the transverse spin compo¬ 
nent. We normalize transverse spin by the magnitude 
of a fully coherent sample with no spin inhomogeneities, 
c{z,t) = S±{z,t)/ and the total coherence in the 
sample is Ctot(^) = f c{z,t) dz. 

Figure [2] shows the typical time evolution of the two- 
domain spin structure. The longitudinal spin component 
exhibits spin domain oscillations as well as diffusion of 
the longitudinal spin gradient, which damps the oscilla¬ 
tions. These measurements reveal that longitudinal do¬ 
mains persist for much longer than expected for classical 
longitudinal diffusion times (^ 25 ms 0) and oscillate 
much slower than trap oscillations (Fig. Efb)). A co¬ 
herent domain wall has a stabilizing effect, which agrees 
with theoretical predictions and previous experiments in 
spin-polarized gases [H, [22l - [2^ . We study this effect in 
more detail by controllably tuning the amount of initial 
coherence in the domain wall. 
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FIG. 2. Time evolution of spin domains with Qnit = 0.74. 
(a) Typical initial spin distribution of |1) (solid line) and |2) 
(dashed line), used to infer S\\. (b) Evolution of S\\. (c) 

Evolution of transverse spin amplitude, c(z,t). Goherence is 
initially confined to the domain wall, in a region less than 
one-third the axial full-width half-maximum of the atomic 
distribution, ~ 0.3 Wax, but rapidly spreads toward the edges 
of the cloud, covering a region of the size ~ 0.7 Wax by 200 ms. 
Ramsey fringes at (d) the trap center z = 0 and (e) z = 
0.3 Wax as shown by arrows in (c), with a sinusoidal fit to 
guide the eye. 


The more striking behavior occurs in the time evolu¬ 
tion of the transverse spin component. Eig. [2fc) shows 
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the spatiotemporal evolution of the amplitude of S± (co¬ 
herence) extracted from Ramsey measurements (Fig.[2jd- 
e)). Ramsey fringe contrast rapidly rises within longitu¬ 
dinal spin domains, where the initial fringe amplitude is 
small (Fig.[2{e)), showing the spread of coherence toward 
the cloud edges. 

Spreading of coherence could result from S± diffusion, 
but the total transverse spin magnitude in the gas in¬ 
creases (Fig. m^d)). This increase implies the appear¬ 
ance of coherence near the cloud edges is not due to 
spin diffusion, but rather to conversion of into S± 
as a result of an instability in the longitudinal spin cur¬ 
rent. This effect has been observed in spin-polarized gas 


systems 


22 


and was described as an experimental 


manifestation of the Castaing instability [12|. Although 
the trapped atomic system possesses different experimen¬ 
tal parameters, the physics governing the phenomenon is 
similar If transverse spin is confined in the do¬ 

main wall, it dephases and the S\\ gradient decays via 
ordinary longitudinal diffusion, since the kinetic equa¬ 
tion decouples for longitudinal and transverse spin com¬ 
ponents [m, [ 2 ^. However, if the ISRE is large enough 
(i.e. large Leggett-Rice parameter), the longitudinal spin 
current becomes unstable, coupling between longitudinal 
and transverse spin dynamics occurs, and the S\\ gradi¬ 
ent decays via transverse diffusion across the coherent 
domain wall. Thus transverse-spin-mediation of longitu¬ 
dinal diffusion increases the timescale for longitudinal dif¬ 
fusion to be comparable with slower transverse diffusion 
times, which are slower due to coherent spin interactions 
absent from longitudinal spin dynamics. 

To study the effect of initial coherence in the domain 
wall, we initialize domain walls with different S± and 
observe their relaxation. Coherence is controlled via the 
same off-resonant laser that prepares the spin domains. A 
short light pulse creates a nonuniform differential atomic 
potential, whose inhomogeneity induces rapid dephasing 
of S± in the domain wall. The degree of coherence is con¬ 
trolled by the pulse length (<0.6 ms) and measured via 
Ramsey spectroscopy. This forced dephasing procedure 
can reduce the initial coherence by over 70%. 

Figure [3fa) shows the longitudinal time evolution for 
different initial degrees of coherence. Both the frequency 
of longitudinal spin domain oscillations as well as longitu¬ 
dinal spin diffusion rates decrease as domain wall coher¬ 
ence is increased. Dynamics for a two-domain structure 
are dominated by the dipole mode, which we isolate via 
the dipole moment of the spin distribution, (z5'|| ( 2 ;, t))^, 
where (..•)^ denotes a density-weighted average to com¬ 
pensate for low signal-to-noise ratio at the distribution 
edges. Figure [3fb) shows evolution of spin-dipole mo¬ 
ments for the coherences in Fig. [3l^a). The frequency / 
and damping rate F of longitudinal oscillations are ex¬ 
tracted from these moments. 

We repeat the experiment for different initial amounts 
of coherence, cinit; calculate the time-dependent dipole 



FIG. 3. (a) Time evolution of S\\ for different initial do¬ 

main wall coherences, from top to bottom Qnit = S±/S^^^ = 
0.74, 0.51, 0.28 respectively in the cloud center, (b) Dipole 
moment time evolution calculated from (a), with a represen¬ 
tative decaying sinusoidal fit to (iii) Cinit = 0.28. (c) Damp¬ 
ing rate and (d) oscillation frequency of the dipole moment 
for different Qnit- Error bars correspond to fit uncertainties 
for dipole moment oscillations and Qnit measurements. The 
shaded band is the result of numerical simulations of Eqn. [T] 


moment; and extract / and F from decaying sinusoidal 
fits to these oscillations. Figure n c) and (d) summa¬ 
rize the results. Both F and / decrease as cinit increases, 
showing the stabilizing effect of a coherent domain wall. 
The oscillation frequency is primarily controlled by two 
competing factors: trapping potential oscillations and 
mean-field-induced spin rotation. In the incoherent limit, 
/ approaches the trapping frequency, as the system ap¬ 
proaches a mixture of distinguishable classical ideal gases 
that diffuse according to classical Boltzmann theory. The 
large cinit limit entrains longitudinal diffusion with slower 
transverse diffusion. 

We compare these measurements with transport the¬ 
ory by numerically simulating Eqn. [T] (see [ 2 ^ for details). 
The shaded regions in Eig. [3fc) and (d) represent one- 
sigma confidence bands from Monte Carlo simulations of 
Eqn.m including statistical fluctuations in n, T, and do¬ 
main wall size, as well as a systematic density calibration 
uncertainty. The data agrees well with theoretical predic¬ 
tions. Discrepancies in / at low coherence is likely due to 
challenges in fitting critically damped oscillations where 
the quality factor drops. Overdamping should occur for 
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FIG. 4. (a) Time evolution of S± for different domain wall 
coherences, from top to bottom Qnit = 0.71, 0.51, 0.30 re¬ 
spectively in the cloud center, (b) Time evolution of S± in 
different spatial regions, denoted by 2 : = (0, 0.3, 0.6) x Wax 
[dotted lines in (i)]. The dashed box indicates a change in 
sampling rate to measure both the fast initial rise of coher¬ 
ence and longer relaxation to equilibrium, (c) Time evolution 
of total ensemble coherence for qnit shown in (a). The solid 
curve is a numerical simulation of Eq.[T]for (i), showing good 
qualitative agreement. Quantitative agreement is limited by 
challenges extracting Ramsey fringe amplitudes in the pres¬ 
ence of noise, dephasing, and small signal in cloud wings. 


Gnit <0.2, but reducing coherence to this level without 
altering the longitudinal spin domains is challenging. 

The effect of the spin instability can be seen more 
clearly in the transverse spin dynamics (Fig. |4]). The 
time evolution of S±{z, t) for different initial domain wall 
coherence is shown in Fig. |4l[a). The rise and spread of 
coherence for different qnit shows similar spatial behav¬ 
ior, but S±_ persists longer when there is more coherence 
in the domain wall initially. Figure |4l[b) shows the evolu¬ 
tion of coherence in different regions of the cloud for the 
high qnit preparation (Fig. 11^(1)). 

Figure |4l[c) shows total ensemble coherence, qot(^): 
for different qnit, calculated by summing the coherence 
across the cloud. Coherence rapidly rises followed by a 
gradual decrease, but with different timescales for dif¬ 
ferent qnit- Maximum total coherence across the cloud 
depends on the initial coherence, and the time to reach 
maximum coherence also increases as qnit increases. 
These results highlight the transverse source of enhanced 
lifetimes of longitudinal spin domains with coherent do¬ 
main walls. The presence of coherence in domain walls 
links diffusion timescales and increases longitudinal dif¬ 
fusion times to be as long as transverse diffusion times. 

Our results validate the quantum Boltzmann equation 



FIG. 5. Instability onset seen in the ratio of maximum en¬ 
semble coherence to initial coherence, c^t^/ 4 o^t, from simula¬ 
tions of Eq. [U plotted logarithmically to highlight instability 
threshold. We vary a and (a) T or (b) £ with other parameters 
set to conditions of Fig. [3l[a)(ii). Since /x oc 1/a, increasing a 
diminishes the instability, as collisional damping grows faster 
than exchange scattering. Decreasing T or increasing i sta¬ 
bilizes the spin current as expected, through increasing adia¬ 
batic rotation during domain wall crossing. Grosses represent 
experimental conditions used here. 

used to simulate the domain wall spin dynamics. As the 
region of experimental interest lies between the collision¬ 
less and hydrodynamic regimes (Knudsen number ^ 1.2), 
no analytic approximation is readily available for the in¬ 
stability criterion, such as was done in [ 2 ^. However, 
there are two primary requirements for driving the in¬ 
stability. First and most important is that the ratio of 
exchange scattering to elastic scattering rates must be 
approximately greater than 1 (quantified by the Leggett- 
Rice parameter fi ^ 3 for this experiment), lest colli¬ 
sions damp spin currents before they can grow large. The 
second requirement is that the number of exchange colli¬ 
sions while traversing the domain wall be small; if not, a 
spin crossing the gradient will smoothly rotate from one 
domain to the other, similar to adiabatic following. 

We extend these simulations to a broader parameter 
space than can be achieved in our system to explore fur¬ 
ther the crossover from classical to quantum diffusion, 
by numerically varying T, n, a, and domain wall size, £. 
Figure [5] shows the effects of T, a, and £ on the spin in¬ 
stability. Because of near-critical damping and spatially 
varying density (thus collision rates), the system exhibits 
a smooth transition between diffusive regimes. Further¬ 
more, since /i is independent of n, density has little effect 
on the instability within the range between collisionless 
and hydrodynamic behavior, and confinement plays little 
role so long as ^ ^ Wax- 

In conclusion, our results reveal the stabilizing effect of 
coherence in a domain wall, leading to an increase in lon¬ 
gitudinal domain lifetimes. These effects are explained 
by coupling between longitudinal and transverse spin dy¬ 
namics, which induces transverse-channel-mediated lon¬ 
gitudinal spin diffusion. While not a direct measure of 
anisotropic diffusion, alteration of diffusive timescales by 
adjusting only the longitudinal-transverse spin coupling 
(via the coherence) is strongly suggestive of anisotropic 
diffusion. An accompanying rise in the total coherence 
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indicates a conversion of longitudinal spin into transverse 
spin as a result of an instability in the longitudinal spin 
dynamics. Longer longitudinal domain lifetimes result 
from more coherence in the domain walls, and suppres¬ 
sion of this instability leads to even longer longitudinal 
diffusion times. Our results agree well with numerical 
simulations and emphasize the importance of quantum 
effects in spin transport in trapped ultracold bosonic sys¬ 
tems, even above quantum degeneracy. 
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